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(— I I Renormalization scheme uncertainties in the next-next-to-leading order QCD pre- 

^ ! dictions are discussed. To obtain an estimate of these uncertainties it is proposed to 

S^ . compare predictions in all schemes that do not have unnaturally large expansion coeffi- 

C^ ' cients. A concrete prescription for eliminating the unnatural schemes is given, based on 

the requirement that large cancellations in the expression for the characteristic renor- 
malization scheme invariant should be avoided. As an example the QCD corrections to 
the Bjorken sum rule are considered. The importance of the next-next-to-leading order 
corrections for a proper evaluation of perturbative QCD predictions is emphasized. 
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In the process of evaluation of perturbative QCD predictions one inevitably encoun- 
ters the problem of renormalization scheme (RS) dependence. As is well known, despite 
the fact that final predictions of a quantum field theory for physical quantities should be 
independent of choice of the renormalization scheme, the results obtained with the truncated 
perturbation series may be numerically different in different schemes. Although the differ- 
ence is formally of higher order in the coupling constant, numerically it may be significant 
in the phenomenological analysis. This effect complicates somewhat comparison of the the- 
oretical predictions with the experimental data. There has been a lively discussion on the 
problem of a proper choice of the RS |]I], ^ ^ ^, ^, ^ . 

However, the problem of a proper choice of the RS is only one of the aspects of the 
RS dependence effect. In fact, the RS dependence of perturbative predictions, instead of 
being an annoying complication, may also become a source of an important information. 
The idea is to estimate the strength of the RS dependence, and to use it as a measure of 
reliability of perturbation expansion. Indeed, we expect the perturbative expression to be 
less reliable when the expansion coefficients are large and/or when the coupling constant is 
not very small. These are precisely the conditions under which we expect also a strong RS 
dependence of perturbative predictions. 

A procedure that is frequently used to estimate the strength of the RS dependence 
consists of calculating the variation of predictions when the renormalization scale is changed 
in some reasonable range. Unfortunately, in the next-next-to-leading order (NNLO) approx- 
imation the freedom of choice of RS is parametrized by two independent parameters, so that 
varying the renormalization scale we scan only a small subset of all the available schemes. 
Also, as has been discussed for example in 0, the concept of a "reasonable" renormalization 
scale is not well defined. 

A more general approach was formulated in 0, in which full freedom of choice of 
the RS is taken into account. This approach is based on the observation that in NNLO 
one may select acceptable schemes using the expression for a RS invariant combination of 
the expansion coefficients. In this note the ideas presented in are further developed. 
A precise mathematical formulation of the constraint eliminating the unnatural schemes is 
given, which allows for a quantitative estimate of the strength of RS dependence for any 
physical observable for which the NNLO corrections are known. The QCD corrections to the 
Bjorken sum rule for polarized structure functions are discussed as an example. The results 



obtained for other observables are also briefly summarized IQ. 

Let us consider a NNLO expression for a physical quantity R depending on a single 
energy variable P: 

R^^\P)=aifi)[l + n{fi/P)aifi) + r2ifi/P)a\fi)], (1) 

where /i dentes the renormalization scale, a{fi) = g^{^)/{4iT'^) denotes the running coupling 
constant that satisfies the NNLO renormalization group equation: 

IJ-r- = -b a [1 + cia + C2a , (2) 

and where 

ri(/i/P)=r!°^+Mn(/i/P), (3) 

r2(/i/P) = rf + (ci + 2ri°))6 ln(/i/P) + (6 ln(/i/P))2, (4) 

The coefficients r> are independent of yU. The commonly used renormalization group im- 
proved expression for R is obtained by assuming that fi = kP, where k is a. constant. In 
the following it is assumed, that the quark mass effects are approximated by a step-function 
energy dependence of the number of active quark flavors. 

The results of perturbative calculations are usually expressed in the modified minimal 



subtraction (MS) scheme 0. Other choices of RS are possible, corresponding to different 
choices of the finite parts of renormalization constants. Results in various schemes are related 
by a finite renormalization, which in our approximation is equivalent to a redefinition of the 
coupling constant: 

(^Wsif^) = a{fi)[l + Ai a{fi) + v42 a^(/i)]. (5) 

The constants Ai and A2 may in principle be arbitrary. The expansion coefficients ri, r2 
and C2 depend on the choice of RS. Also the parameter A does depend on the RS [|l^]. The 
relevant formulas describing this dependence are collected for example in 0. However, in 
NNLO there exists an RS independent combination of the expansion coefficients |^, [Tl|, p!2| : 



P2 = C2 + r2 -ciri -r^ (6) 

As was pointed out in [^ , a convenient description of the RS dependence of P*^^^ is obtained 
when ri and C2 are chosen as independent parameters characterizing the choice of RS. The 



parameter ri controls the degree of freedom in the approximants which is equivalent to the 
freedom of choice of the renormalization scale /i, and C2 controls - in the terminology of 
Stevenson [Q] - the dependence on the choice of the renormalization convention. With this 
parametrization the NNLO RG- improved expression for R may be written in the form 0: 

i?(2)(p/A^;r-i,c2) = a(P)[l + ria(P) + r2(ri,C2)a2(P)], (7) 

where r2 is expressed as a function of ri and C2 using the Eq. ^. The numerical value of the 
running coupling constant is obtained from the implicit equation: 

bH^)=rY^-n + ^^'\a,C2), (8) 

where 

$(2)(a,c2) = ciln(;— ) + - + Ciln(cia) + 0(a), (9) 

2ci a 

which is obtained by integrating the Eq. ^ with an appropriate boundary condition. The 
explicit form of $(a, C2) may be found for example in |13|. The presence of Ap^ in the 



Eq. I is a result of taking explicitly into account the one loop Celmaster-Gonsalves relation 



between A parameters in different schemes |T0[ , which is valid to all orders of perturbation 
expansion. 

The fact, that numerical value of the predictions depends on the choice of RS, stimu- 
lated the interest in prescriptions, which pick up schemes distinguished by some arguments 
of physical or mathematical character. Several such prescriptions for selecting a "good" RS 
have been discussed, like the prinicple of Fastest Apparent Convergence (FAC) [|I|], the Prin- 
ciple of Minimal Sensitivity (PMS) 0, the momentum subtraction prescription |^, ^, and 
other m. However, regardless of our choice of the optimal scheme, there is always a contin- 
uum of schemes which are close to the one preferred by us. Since any choice of a preferred 
RS may have only a heuristic motivation, the predictions obtained in schemes close to the 
optimal scheme also must be somehow taken into account. A natural way to accomplish this 
is to supplement the prediction in a preferred scheme by the estimate of the strength of the 
scheme dependence. 

In order to obtain a meaningful estimate of the strength of RS dependence one should 
calculate variation of predictions over all a priori acceptable schemes. To this end we must 
provide some well motivated criteria for selecting the admissible scheme. The use of some 



selection criteria cannot be avoided, because by taking variation of predictions over all pos- 
sible schemes we would obtain a meaningless result. The reason for this is obvious - clearly 
one may make a bad choice of scheme, in which the perturbation expansion could have 
unnaturally large expansion coefficients, and one may artificially obtain very large higher 
order corrections and very strong RS dependence, without any physical significance. As was 
pointed out in [0 , in order to eliminate unnatural schemes one may use the expression for the 
RS invariant p2- This invariant appears as an NNLO expansion coefficient in the manifestly 



RS-invariant evolution equation for a physical quantity [0, |T2|, |14[ , and therefore it may be 
taken as a natural RS-independent characterization of the magnitude of the NNLO correc- 
tions. The observation made in 0] was that one should avoid using schemes that introduce 
large cancellations in the expression for p2- 

To make this condition more concrete let us introduce a function which describes the 
degree of cancellation in the expression for p2 in a quantitative way: 

cr2(n,^2,C2) = |C2| + IrsI +Ci|ri| + r^. (10) 

The requirement, that the contributions to p2 should not be too large compared to the 
magnitude of p2 itself, may be now written in the form: 

cr2{ri,r2,C2) < l\p2\, (11) 

where / > 1 is some constant . The constant / determines how strong cancellations in the 
expression for p2 we want to allow. For arbitrary / it is straightforward to obtain an explicit 
description of the region of allowed parameters in the (ri,C2) plane. In the case p2 > cf/A, 
which is what will be needed for the Bjorken sum rule, let us define: 
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Vp2(/-1)/2, (12) 



r™- = [_Ci + v/cf + 2(/-l)p2]/2, (13) 

cr" = -P2{1 - l)/2, (14) 

c-- = p2(/ + l)/2, (15) 

4'^* = cir™" + c^"^\ (16) 

For C2 > the region of allowed parameters is bounded from above by the line joining the 
points (r^™, 0), (r™", 4"*), (0, c^*"^), {r"^"'', c^"^), (r™'^^, 0). For C2 < the region of allowed 



parameters is bounded from below by the lines: 

C2(ri) = rl + c™" for r™" < ri < 0, (17) 

C2(ri) = r2 + ciri + cr" for < n < r^''. (18) 

It is interesting that the FAC scheme, corresponding to ri = 0, C2 = P2, always belongs to 
the allowed region. We also see, that for P2 = - i.e. null NNLO correction - the only 
scheme satisfying Eq. |ll] is ri = 0, C2 = (and of course r2 = 0). 

For a final specification of the allowed region in the (ri, C2) plane we have to choose 



some value of / in the Eq. |TT]. To make a meaningful choice let us observe, that it is a 
rather natural requirement to have the PMS scheme in the allowed region. The PMS scheme 
parameters are determined by the condition: 



An approximate solution of these quations has the form [15 



r^^s ^ o(a^*^^), cf^^ = ^p2 + 0(a^^^). (20) 

It is easy to see that taking / = 2 we would have the PMS scheme right at the boundary of 
the allowed region in the (^1,02) plane. 

Therefore, in order to have an estimate of the strength of the scheme dependence in 
NNLO we propose to calculate the variation of predictions using Eq. |^ and Eq. |], when the 
scheme parameters ri and C2 are changed within the region satisfying the Eq. |ll] with / = 2. 

Let us now consider the QCD correction Rbj to the Bjorken sum rule for the polarized 



structure functions in the deep-inelastic scattering [16 



/' dx[gl{x,Q')-gl{x,Q')] = \i^[l~RBAQ')]- (21) 

Jo D Qy 

This sum rule has been recently a subject of intense experimental [|T^ [1^ and theoretical 



1^, |20| investigation. The QCD perturbation expansion for Rbj has the form of Eq. |l], with 

P^ = Q^ ■, where the expansion coefficients in the MS scheme with four quark flavors are 
given by [|2l|, |T9|: 

rf^ = 3.25 rf^= 13.8503. (22) 



Together with ci = 1.54, cf ^ = 3.0476 this imphes p2 = 1.3304 (6 = 25/6). In Figure 1 we 
show R^j as a function of ri and C2 for P/Ajjg = 8.6 (this corresponds to < Q^ >= 5GeV'^ 
- as used by the SMC Collaboration to present the combined experimental data |T^ - and 



,(4) 



^J¥s ~ 0-26 GeV as preferred by the Particle Data Group [^), together with the region 

(2) 

of allowed values for the scheme parameters. We see that R^j has indeed a saddle point, 
representing prediction in the PMS scheme, close to the / = 2 allowed region. We have: 

i?£f ^(8.6) = i?g](8.6; 0.0852, 2.1302) = 0.126525. (23) 

The estimate of the strength of the scheme dependence is obtained by comparing the pre- 
dictions inside the allowed region. We find that the extremal values are attained at the 
boundary of the allowed region: 

i?^J^(8.6) = R%{S.6; -0.668, -0.219) = 0.12697, (24) 

i?™"(8.6) = i?g](8.6; 0.352, 0) = 0.12597. (25) 

For comparison we also show in Figure 1 the range of scheme parameters satisfying the 
constraint given by Eq. ^ with / = 3, in which case a small neighbourhood of the PMS 
parameters is taken into account. In this case we obtain: 

i?™J^(8.6) = R%{S.6; -0.838, -0.628) = 0.12722, (26) 

i?™"(8.6) = i?g](8.6; 0.617, 0) = 0.12551. (27) 

In Figure 2 we show the prediction for i?^ • as a function of P/Ajjg near the point 
P/Ajj^ = 8.6. This figure shows how the uncertainty of the predictions arising from the 
RS dependence may affect the fit of Aj^ to possible experimental results. We see that the 
obtained scheme dependence appears to be very small. In fact it is much smaller than the 
uncertainty arising from the estimated error in the fitted value of A-jjg - the Particle Data 

Group gives Aii^ = 0.316 — 0.214:VeV, which for Q"^ = 5 GeV^ gives variation of Q'^/A'^ in 

the range (7.08)^ — (10.45)^. The scheme dependence uncertainty is also much smaller than 
the experimental errors in the presently available experimental data [|l^, |T^. This shows 
that the nj = 4 NNLO formula for the QCD corrections to the polarized Bjorken sum rule 
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seems to be a potentially good source of precise information on A;p^ once the experimental 
accuracy is improved. 



Let us also note that the MS coefficients lie outside the allowed region indicated 



in Figure 1, i.e. in our approach the MS scheme has to be considered as a scheme with 
unnaturally large expansion coefficients. Also, Rj^f{8.6) = 0.1f969 lies outside the range of 
values obtained in our estimate of scheme dependence uncertainties. However, the difference 
between R^f and Rsf^ is still rather small for Q'^/Aj^ close to (8.6)^. 

For the sake of comparison it is interesting to consider the RS dependence of the 
NNLO expression for Rbj for Uf = 3. Taking three quark flavors we have: 

rf^ = 3.5833 rf^ = 20.2153, (28) 

which together with ci = 16/9 and c^'^ = 4.471 implies p2 = 5.4757 {b = 9/2). This is 
a significantly larger value than for four flavors. The dependence of i?^ • on ri and C2 for 
Uf = 3 for P/Aj^ = 8.6 is shown in Figure 3. We flnd: 

rP^s ^ i?g](8.6; 0.33, 8.83) = 0.121634. (29) 

Repeating the procedure used for n/ = 4 we obtain for the / = 2 allowed region: 

i?™«^(8.6) = i?g](8.6; -1.655, 0) = 0.12748, (30) 

i?™" = i?g] (0.989, 0) = 0.11757. (31) 

For the / = 3 allowed region we flnd: 

/2™«^(8.6) = i?g](8.6; -2.340,0) = 0.13412, (32) 

i?™"(8.6) = R%{S.6; -2.340,6.791) = 0.10683. (33) 



Note that also in this case the MS coefficients lie outside the allowed range for scheme 
parameters. We have R^f = 0.11102, which is close to the lower limit of values obtained 
by varying over the / = 3 allowed region. We see that the theoretical uncertainty in the 
Hf = 3 formula is larger than in the nj = 4 formula. This is a result of stronger dependence 

of i?*^^-* on the scheme parameters in this case, and of the larger range of parameters which 



may be considered to be natural. Using our method we may now quantitatively compare 
reliability of the Uf = 3 and Uf = 4 predictions for Rbj ■ Concerning the energy dependence 
of the theoretical uncertainties in this case let us note, that for lower values of Q'^/A'^ the 
RS dependence of the three-flavor expression, estimated according to our prescription, is 

(2) 

considerably larger than that found above. The problem of the low energy behavior of R^, 
and a possible method for reducing scheme dependence at low energies would be discussed 



in a separate note 23 



We conclude this note with several remarks. 

1. It has to be emphasized that our estimate of the strength of the RS dependence 
does not provide an error estimate in the mathematical sense, i.e. there is no theorem that 
would guarantee that the true result lies within the obtained range of variation for the 
prediction. The method described above provides only a consistent framework for making an 
unbiased educated guess about possible theoretical precision of the prediction. The method 
avoids unnecessary constraints on the considered renormalization schemes, it eliminates from 
analysis the unnatural schemes, and it is sensitive to the magnitude of the NNLO correction 
via dependence on p2. The estimate of the scheme dependence uncertainty obtained in this 
way should be particularly useful for a quantitative comparison of reliability of perturbative 
predictions for different physical quantities, evaluated at different energies. In fact one 
may reasonably expect that thus obtained relative estimate of the reliability of perturbative 
predictions would be insensitive to certain freedom in choosing I in the Eq. |ll|and to possible 
variations in the choice of the function o"2. The proposed method should be also useful in 
determining the regions of applicability of perturbation expansion. Indeed, since we have 
taken great care to eliminate schemes with unnaturally large expansion coefficients, in order 
to avoid introducing any spurious RS dependence, then any large variation of predictions 
obtained according to our method would be an unambiguous sign of a limited applicabihty 
of the NNLO expression in the considered energy range. 

2. The RS invariant p2 plays a fundamental role in our approach. However, Stevenson 
in his anlysis of RS dependence used an expression for the NNLO invariant that differs 
from the one used here by a constant. We think that the expression used in this note is 
more natural for characterization of the magnitude of NNLO corrections because it appears 



in the RS invariant evolution equation for a physical quantity [|TT|, [T2|, [l^. To emphasize 



the RS independent character of the evolution equation obtained in [ITl], |l2l ^M let us recall 



the argument indicated in [^. Let us take i?'-^^ as given by Eq. |I] in arbitrary scheme, 
differentiate over P, then express a{P) in terms of R^'^\ and finally express P(dR^'^ydP) as 
a series expansion in R^'^\ We obtain: 

jp{2) oo 

P^ = -6(i?^'))'[l + cii?(^) + P2(i?(^))^ + E MR^'^n (34) 

fc=3 

The first three expansion coefficients in this equation are RS independent and coincide 



with the coefficients obtained in |11, O, 14 1. Higher order expansion coefficients in this 



equation are scheme dependent, which is a result of the fact that we have only used NNLO 
approximation for a physical quantity. (If higher order approximant would be taken, then 
more coefficients in the above equation would become independent of the choice of RS.) This 



shows, that the equation obtained in [|ll| [I2|, |14| is indeed RS independent, and that it is 



not a result of a particular choice of the scheme. 

Another argument supporting our choice is the fact that p2 naturally appears in the 
leading order solution of the PMS equations (Eq. 



3. The fact that we base our method on the existence of the RS invariant p2 has an 
obvious consequence, that it cannot be applied to quantities for which only next-to-leading 
(NLO) corrections are known. This should not be surprising however, because the NLO 
approximation has a rather special status in QCD. In fact, the NLO correction allows only 
to set the proper energy scale for the prediction. For example, using the FAC scheme we 
obtain in NLO an identical energy dependence for all possible physical observables depending 
on one energy variable, except for a shift on a logarithmic energy scale. It is only at NNLO 
that physical quantities receive process-specific corrections. It is therefore natural that only 
at NNLO and beyond one may test reliability of the perturbation expansion. This indicates 
a fundamental importance of the NNLO calculations for a proper comparison of theoretical 
predictions with experimental data. Let us also note that it is straightforward to extend our 
approach to the case of four-loop or even higher order corrections if such calculations would 
ever be done. 

4. Let us discuss some other methods of estimating the reliability of perturbative 
predictions. A simplest alternative is to compare NLO and NNLO predictions, taking the 
difference as an estimate of the theoretical uncertainty. Unfortunately, it is obvious that such 
an estimate would be strongly dependent on the choice of scheme. One may also try to guess 
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the magnitude of higher order corrections. Again, the problem here is the RS dependence 
of such an approach. Note also that the asymptotic expression for higher order corrections 
obtained in P3, which has been sometimes used to evaluate the theoretical uncertainty, is 



known to be erroneous pq]. Finally, in 27 it is proposed to obtain theoretical uncertainty 



by comparing predictions obtained in the PMS, FAC and MS schemes. Unfortunately, 
this method also may give misleading results, since the predictions in the PMS and FAC 
approach are known to lie rather close to each other regardless of the magnitude of the 



NNLO correction, and the choice of MS in such procedure is completely arbitrary, without 
any theoretical or phenomenological motivation. In fact, as was shown above, comparing 



PMS and MS we may obtain larger uncertainty than that suggested by our method. 

5. The method described in this note has been applied to other quantities for which 



the NNLO expansion is known pSl E3, ES 



The results for the Gross-Llewellyn-Smith sum rule are similar to those obtained 
above for the Bjorken sum rule. The NNLO expression for nj = A appears to have small 
scheme uncertainty, while the Nf = 3 expression seems to be less precise, particularly at 
lower energies. 

In the case of the QCD corrections to the tau lepton decay we confirm the observa- 
tions made in [0 . By varying predictions over the set of schemes allowed by the constraint 
given by Eq. [Tl| we find that the frequently used NNLO perturbative expression is unstable 
against the change of RS. (In the condition on the allowed schemes was not specified 
in a quantitative way.) However, if one takes a contour integral representation for i?^, one 



obtains a stable result |2^, with theoretical uncertainty of the same order as the present ex- 
perimental uncertainty. Neglecting nonperturbative contributions for simplicity and fitting 
^M5 ^° ^^^ experimental central value R^^^ = 0.2 we obtain 0.396 Gel^ in the generalized 
PMS method, and 0.386 — 0A15GeV variation for the / = 2 allowed region for the scheme 
parameters. 

Application of our method to the case of the Uf = 5 QCD corrections to the e"^e~ 
annihilation into hadrons reveals a surprisingly large RS dependence, despite rather high 
characteristic energy scale. The large scheme uncertainty is a consequence of a large value 
of the RS invariant p2- However, a closer analysis shows that large value of this invariant is 
due mainly to the effect of analytic continuation from the spacelike to timelike momenta. A 
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proper treatment of the analytic continuation results in a considerable improvement of the 
stability of the predictions with respect to change of the RS |25]. Similar remarks apply to 



the case of the QCD corrections to the Z-boson decay rate to hadrons, which is however 
more complicated because of the presence of mass dependent corrections. 

Summarizing, we may say, that a method has been proposed for analyzing the un- 
certainties in the NNLO predictions, introduced by arbitrariness in the choice of the renor- 
malization scheme. This method is based on a condition that eliminates schemes which 
have expansion coefficients giving unnaturally large cancellations in the expression for the 
characteristic NNLO RS invariant. Application of the method has been illustrated using as 
an example the QCD corrections to the Bjorken sum rule for polarized structure functions. 
It was found that for four quark flavors these corrections have small scheme dependence 
uncertainty. Our considerations show the fundamental importance of the NNLO corrections 
for a proper evaluation of the perturbative QCD predictions. 

Author is grateful to the High Energy Physics group at the Argonne National Labo- 
ratory for hospitality and discussion. 
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Figure Captions 



Fig. 1. Contour plot of the NNLO prediction for Rbj as a function of ri and C2, for 
Uf = 4 and Q^/Ajr^ = (8.6)^. The region of values of the scheme parameters satisfying the 
constraint given by Eq. |ll| is indicated for / = 2 (smaller region) and for / = 3. The dashed 
line corresponds to R^j = 0.126526. 

Fig. 2. The NNLO prediction for Rbj as a function of P/Ams (^^ = <5^)- The thick broken 
line in the middle denotes the PMS prediction. The thin solid lines denote the variation of 
the predictions when ri and C2 are changed within a region defined by Eq. |ll| with / = 2. 
Thin broken lines correspond to / = 3. 

(2) 
Fig. 3. Same as in Fig. 1 but for Uf = 3. The dashed line corresponds to R^j = 0.121635. 
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